Anisotropic fluctuations in turbulent sheared flows 
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An experimental analysis of small-scales anisotropic turbulent fluctuations has been performed 
in two different flows. We analyzed anisotropic properties of an homogeneous shear flows and of a 
turbulent boundary layer by means of two cross-wire probes to obtain multi-point multi-component 
measurements. Data are analyzed at changing inter-probe separation without the use of Taylor 
hypothesis. The results are consistent with the "exponent-only" scenario for universality, i.e. all 
experimental data can be fit by fixing the same set of anisotropic scaling exponents at changing 
only prefactors, for different shear intensities and boundary conditions. 

I. INTRODUCTION 

Statistical theory of turbulence is often focused on homogeneous and isotropic flows [1]. Experimentally, however, 
we know that isotropy holds approximately, with different degrees of justification, depending on the geometry of the 
boundaries and on the driving mechanism. Therefore, a realistic description of turbulence cannot ignore anisotropic 
and non-homogeneous effects, especially in regions close to the boundaries, and/or at scales close to the integral scale, 
Lq, where the injection mechanism can strongly affect velocity fluctuations. The interest of quantifying anisotropic 
and non-homogeneous effects is also linked to the important issue of "recovery of isotropy", i.e. the problem of 
"small-scales universality" . Surprisingly enough, recent experimental and numerical works [2-7] have detected the 
survival of anisotropic turbulent fluctuations till to the Kolmogorov scale, r/. These findings have stimulated a lot 
of further experimental, numerical and theoretical work focused to develop proper analytical tools [8] and to extend 
the available experimental/numerical data sets [2-7, 9-12]. Many progresses have been done. For example, the so- 
called puzzle of "persistence of anisotropics" in small-scales -high Reynolds numbers- sheared flows, has been recently 
understood as the effect of the presence of strong anomalous anisotropic fluctuations [7, 13]. The attention is mainly 
focused on correlation functions based on gradients (to probe Reynolds number dependencies) or to the projections on 
isotropic/anisotropic sectors of multi-points velocity correlation functions, S^ a '(ri, . . . ,r„) = (v ai (ri) . . .v an (r n )), 
where we use {a} as a shorthand notation for the ensemble of indexes a\, . . . ,a n . When all spatial separation, 
r\ . . . r„, are in the inertial range, rj <C |t"^ — Tj\ ~ r <C Lq, one expects the existence of power laws behavior 
under a uniform space dilation: S^(Xri, . . . , Xr n ) = A""J (t*i, . . . , r n ). Most of the recent works in anisotropic 
turbulence concentrated on determining the values of the exponent, £(n), as a function of the order of the correlation 
function, n, and of its anisotropic properties. Indeed, an important step forward has been done by realizing that 
different projections of the multi-point correlation functions on different irreducible representation of the group of 
rotation, SO (3), possess different scaling properties. The idea is to decompose any correlation function in a complete 
basis of eigenfunctions with defined properties under rotations. Each eigenfunction identifies a specific anisotropic 
sector with total angular momentum, j, and its projection, m, on a given axis. It is believed that the scaling properties 
of the projections on different sectors possess different scaling exponents, £,j{n). 

Exponents for the fully isotropic sectors are labeled by j = while more and more anisotropic fluctuations are 
measured by higher and higher values of j [8, 10]. 

The higher-than-expected presence of small-scales anisotropic fluctuations raises questions about their universal or 
non-universal origins. In other words, one is interested to control if all flows posses the same, or similar, anisotropic 
small-scales fluctuations independently on their large-scale behavior. Of course, full universality cannot be expected, 
one is tempted to believed to a "exponent-only" scenario, i.e. only the scaling exponents, S,j{n), pertaining to each 
different anisotropic sector, are universal, while the overall correlation functions intensities are not. This hypothesis 
is inspired by both theoretical reasons [8, 14] and similarities with what observed for isotropic fluctuations [15]. 

Up to now, there are a few experimental and numerical data sets where universality of anisotropic fluctuations 
has been probed. In particular, as of today, we have some detailed experimental investigation of anisotropic small- 
scale fluctuations in homogeneous shear flows [3, 4], in atmospheric boundary layer [10-12] and in windtunnel data 
[16]. From the numerical side, only a few Direct Numerical Simulations (DNS) in highly anisotropic flows have been 
performed with the aim to explicitly test the small-scales properties of anisotropic fluctuations [6, 9, 13, 17, 18]. 
The situation is still moot. On the experimental side, because of the difficulty to have multi-point multi-components 
velocity measurements one can access only the j = 2 sector. On the other hand, DNS can properly disentangle 
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fluctuations of all sectors, but due to limitations in the Reynolds numbers, only results in the j = 4 and j = 6 sectors 
have been obtained with some accuracy. The j — 2 sector in the numerical works [6, f 7] was not measurable due to 
strong finite Reynolds effects. Results from different experiments, with different geometries and different large scale 
structures, are in fairly good agreement concerning the j — 2 sector up to moment n — 6. Putting together all results 
of numerics and experiments one recovers a scenario for anisotropic fluctuations consistent (not in contradiction) with 
the "exponent-only" picture of universality. Still, more tests in both experiments and numerics are needed. 

The aim of this paper is to present a new systematic assessment of anisotropic fluctuations in sheared flows at 
changing both the experimental set-up and the shear intensity, fn particular, we have measured small-scales turbulent 
properties in a homogeneous shear flow (HS) and in a turbulent boundary layer (TBL). One of the novelties here 
presented consists in the using of measurements from two cross-wire probes at changing their separation, i.e. we do 
not need to use Taylor hypothesis of frozen turbulence in order to extract information at different scales. This kind of 
multi-points multi-components measurements are necessary in order to disentangle contributions from isotropic and 
anisotropic fluctuations and among different kind of anisotropic fluctuations. 

Our results support the "exponent-only" scenario. We found good qualitative, and quantitative, agreement of the 
anisotropic scaling exponents in both HS and TBL flows. Moreover our results are in agreement with the previously 
measured values in different experiments with different Reynolds numbers and different shear intensities. 
The paper is organized as follows, fn sec. II we present the details of the two experimental apparatus including some 
typical large-scale measurements to validate the laboratory set-up. In sec. Ill we present the scaling properties for 
both HS and TBL flows. Conclusions follow in sec. IV. 



II. EXPERIMENTAL SET-UP 

The data we are going to discuss concern two different experiments, both conducted in the 1.30 m x 0.90 m test 
section of a 7 m long open return wind tunnel. The first data-set has been obtained in a nominally homogeneous 
shear flow (HS), characterized by a constant velocity gradient. The second one refers to measurements performed in 
the logarithmic region of a zero-pressure gradient turbulent boundary layer (TBL). 

A. Homogeneous shear flow 
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FIG. 1: Homogeneous shear flow: mean velocity profile as a function of the non-dimensional coordinate in the mean gradient 
direction, y/h. In the inset, the normalized streamwise turbulent intensity u rlnB . Velocities are normalized in terms of the 
centerline velocity, U c = 10.2 ms _1 , while the test section height is h = 0.9 m. 

The set-up of the homogeneous shear flow is based on the original idea proposed in [19]. The mean shear is generated 
with a device consisting of a series of 15 adjacent small channels, equipped with screens of different solidity to produce 
suitable pressure drops. The channels are followed by flow straighteners in such a way that sufficiently downstream 
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TABLE I: Basic parameters for the homogeneous shear flow. Symbols are defined as follows: The rms fluctuation intensity in 
the streamwise direction is u rms = ((it — U) 2 ) 1/>2 , and analogous definitions apply for v lnls and m lms . p uv = (wu)/(w rms v Ilna ) 
is the correlation coefficient, e is the energy dissipation rate evaluated in terms of the one-dimensional spectrum Eh(ki): 
15 v J kiEh(ki) dni . £t is the transverse integral length scale, J ( u(x, y, z) u(x, y,z + r z ) ) dr z /«r ms and r\ = (u A /e) 1 ^ 4 is the 
Kolmogorov scale. The Taylor- Reynolds number is Rex = Aw rms /V, where the Taylor scale follows from ^/15^u? ms /e and the 
shear parameter is S* = \J Su? ms /e. 

the mean velocity profile U(y) is linear in the core region. The data shown in figure 1 correspond to a measurement 
approximately 4.8 m downstream of the apparatus, where the flow is already well developed. Concerning fluctuations, 
the deviations from the ideal constant profile of the streamwise turbulence intensity u rms {y) = ((u — U) 2 )? are 
comparable with what observed in similar set-up [2]. In particular, in the central part of the test section where the 
data discussed below have been acquired, they are of the order of 7 %. The dimensionless shear rate S* ~ 5 (table I), 
is a factor two smaller than what achieved in the logarithmic part of the turbulent boundary layer (table II). 

B. Boundary layer 

The boundary layer develops on the smooth surface of the lower wall of the tunnel, where a nominal zero pressure- 
gradient is achieved by adjusting the upper wall. The measurements have been performed on the centerline of the 
test section, 6.0 m downstream the tripping device at the end of the contraction. With an external velocity Uoo of 
11.5 m/s, the thickness of the boundary layer at this location is approximately S ~ 40 mm, while the Reynolds number 
based on the momentum thickness Reg is approximately 6500, well within the range pertaining to a fully developed 
turbulent boundary layer. The friction velocity, u T — \Jr w j p with t w the average shear stress at the wall and p the 
constant density, estimated from the mean velocity profile with a Clauser chart, is found to be u T — 0.43 ms -1 , in 
good agreement (±8%) with a direct measurement by means of a Preston tube. The streamwise mean velocity U(y) 
and the fluctuation intensity u rms (y) profiles are displayed in Figure 2. Both curves show that the flow complies to 
the requirements of a fully developed turbulent boundary layer. 
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FIG. 2: Turbulent boundary layer: Mean velocity profile U + vs wall normal distance y + in wall units. The dotted line 
corresponds to the linear behavior U + = y + , while the solid line indicates the log-law: U + = (l/ft)lnj/ + + B, with k = 0.41 
and B — 5.1. The inset shows the streamwise turbulent fluctuation u^ ms . Wall units are defined in terms of the friction velocity 
and the kinematic viscosity u, thus U + = U/u T , uf ma = u TlnB /u T and y + — u T y/v. 
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TABLE II: Basic parameters for the boundary layer. For definitions see caption of table I. 



C. Data acquisition 

The instrumentation consists essentially of a couple of sub-miniature X-probes, mounted on streamlined supports 
in order to minimize interference effects, and separated in the transverse direction, see figure 3. The wires are 2.5 fi m 
in diameter, 0.5 mm in effective length and in separation, oriented at ±45° with respect to the streamwise direction. 
They are operated at an overheat ratio 1.9. Single component sub-miniature probes (diameter 2.5 ^m, length to 
diameter ratio equal to 200 to minimize conduction losses) have also been used to measure the velocity profiles. 
The signals from the two X-wires are simultaneously digitized at 21 kHz with a 16-bit data-acquisition board, after 
being low-pass filtered at the Nyquist frequency and suitably amplified to achieve a good signal-to-noise ratio. The 
frequency response of the hot-wires, measured in the free stream at a reference velocity U = 6 ms _1 , is larger than 
15 kHz. In this way, the resolution needed to analyze the small scales behavior down to the Kolmogorov length-scale 
i] - typically of the order of 0.1 -j- 0.2 mm - is guaranteed. The X-wires are calibrated in situ against a Pitot tube. 
The mapping between the output voltages and the components of the velocity (u,v) is obtained by varying both the 
reference velocities and the orientation of the probe with respect to the flow (see, e.g., [20]). The calibration was 
repeated at the end of each set of measurements, to check that no voltage drift had occurred. 

As for the length of the signals, they consist typically of 4 -f- 8 x 10 6 samples, corresponding roughly to 2 x 10 4 eddy 
turnover times for both flows. Convergence of the statistics has been checked up to the 6 th moment. 

III. RESULTS 

Scaling properties of anisotropic fluctuations are traditionally addressed through objects that are identically zero 
in homogeneous isotropic conditions. Typically, the study has been confined to the co-spectra [21]. Recently a new 
extended set of observables has been proposed in the context of the SO (3) decomposition. The idea is to exploit the 
expansion of any generic statistical observable in terms of a suitable eigen-basis with a well characterized behavior 
under rotations [4, 11]. Let us focus, for instance, on the generic element in the space of second order correlation 
tensors 

S<xi<x2 ( r ) = ( Sv ai (r)5v a2 (r) ) , (1) 

where Sv a (r) denotes the a th component of the velocity increment at two points separated by the vector r, Sv a (r) = 
v a (x + r) — v a (x). The appropriate SO(3) decomposition of (1) reads [8]: 

oo +j p(j) 
3=0 m=-j q=l 

Here the index j denotes a sector, to be understood as a subspace invariant with respect to rotations, denote 
the appropriate basis function [8] which depend on the unit vector f and p(j) counts the number of irreducible 
representations. In particular, j — labels the isotropic sector, while sectors of increasing anisotropy correspond to 

(2) 

higher and higher j's. Information on the dynamics of the system is now captured by the coefficients Sj^ which 
depend only on distance r. 

The invariance under rotations of the inertial terms of the Navier-Stokes equations suggests that small-scales 
statistics depend only on the sector under consideration. For Reynolds large enough, scaling laws of the projection 
are therefore expected in the form 
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n separation Observable 

2 n Si 2 (n) 

4 ri Sm 2 (ri) 5'i222(ri) 

6 ri S , iiiii 2 (ri) 5*111222 (ri) £122222 (ri) 

2 r 3 Si 2 (r 3 ) 

4 r 3 Siii 2 (r 3 ) S'i222(r' 3 ) 

6 r 3 S , imi2(r 3 ) 5*111222 (r 3 ) S'i 22222 (r 3 ) 



TABLE III: List of observables with null contribution from the j = and 1 sectors. 

where the scaling exponent explicitly depend on the sector, j, while the argument 2 reminds that we are presently 
dealing with a second order tensor. The machinery can be easily extended to structure functions of any order n, 

S ai ... an (T) = (Sv ai (r)...Sv an (v)), (2) 

whose projection on the proper SO (3) basis will possess a scaling behavior: 

In this context, the recovery of isotropy at smaller and smaller scales correspond to the existence of a hierarchy of 
exponents £j = o( n ) < £j( n ) [7, 8]. 

Scaling laws for the anisotropic sectors has been recently addressed by using different DNS databases [6, 7, 18]. 
From the experimental point of view, the evaluation of the proper SO (3) components of a given correlation tensor is 
hampered by the limited information on its spatial dependence. In the latter case, the simplest approach is to make a 
selection of tensorial components such as to cancel out the isotropic contribution in the expansion (2). For example, 
the component <Si2(r) in the direction r — (0,0, r 3 ) vanishes in a purely isotropic ensemble. Still, in principle all 
anisotropic sectors may influence its behavior. One may follow two possibilities. Either one may extract the whole 
anisotropic spectrum by making a multi-parameter fit in all sectors [10] or may assume, as done in the present paper, 
that at scales small enough the correlation function is dominated by the leading anisotropic contribution [4, 11]. 
Considering that in the geometrical set-up of our interest the j = 1 sector is absent by symmetry, one assumes that 
in the small-scales limit (at high Reynolds numbers) the leading behavior of (2) is given by the j = 2 sector: 

+ 2 p(2) 

s aia2 (r) -EE 

m=-2 9=1 

Clearly, using this procedure, systematic, out-of-control, errors are introduced by neglecting the higher j sectors. 
Similar consideration can be extended to tensorial correlation functions of any order. For example, table III 
lists several observables which, according to the previous discussion and the symmetries of the experimental 
set-up sketched in figure 3, do not present contributions from both sectors j = and j = 1. In the table, the 
suffixes 1, 2 and 3 correspond to direction x\ = x, x^ = y and x$ = z, respectively. The objects reported in the 
first three lines depend on the streamwise separation n, they can be evaluated by using a single X-wirc probe 
and Taylor hypothesis. Those on the last three lines depend on the transverse separation r 3 and can be mea- 
sured only by using at least two X-wire probes. Hereafter we mainly present results based on the two-probes approach. 

Considering the schematic of figure 3, the two points measurements consist in the acquisition of u and v at two 
points separated in direction z (see caption). 

On the other hand, for example, the single point measurement of Sm2(ri) with u\ = u and U2 — v, as a function 
of time yields 

Sim^i) = 

([u(t + n/U) - u(t)} 3 [v(t + n/U) - v(t)}) . (3) 

This approach has been used e.g. in [4] in the context of the homogeneous shear flow and in [11] at a single location 
in the atmospheric boundary layer to address the scaling properties of the j = 2 sector. In [12], two single component 
wires, at fixed separation in the transverse direction z, are used in connection with Taylor hypothesis to extract the 
scaling exponent of the j = 2 sector while a similar procedure with two X-wires separated in the shear direction y 
permits to investigate the scaling behavior of the j = 1 sector. 

Purpose of the present work is to by-pass the use of Taylor hypothesis by using the configuration described in the 
schematic of figure 3. This allows to compute the anisotropic observables depending on r 3 (table III) by continuously 
changing the transverse separation between the two probes. 
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FIG. 3: Schematic of the measurement configuration: each of the two X- wires, separated in the transverse direction, z = £3, 
detects two velocity components u and v (in the streamwise, x = xi, and in the shear direction y = X2, respectively). The 
separation between the two probes was measured by means of a CCD camera with a high magnification lens and could be 
changed by means of a precision transverse gear from 0.8 mm to 70 mm. 



A. The homogeneous shear flow 



The global parameters characterizing the homogeneous shear flow are summarized in table I. In order to allow for 
the direct comparison of the data for the homogeneous shear flow with those for the boundary layer that are presented 
in subsection III B, a common normalization procedure is used. For the homogeneous shear the relevant characteristic 
velocity is defined as 



«r = (4) 

where the total shear stress is given by t = vS — (uv) where 5 = dU/dy is the mean shear. 

In figure 4, the second order mixed structure function 512(7*3) is plotted as a function of the transverse separation. 
In addition to the r 2 behavior at smaller scales and to the large scale saturation, a power-law at intermediate 
scales emerges distinctly, allowing us to measure the scaling exponent £2(2) with good accuracy. The estimate 
£2(2) = 1.2 ± 0.07 is indicated by the solid line, while the horizontal plateau in the inset, displaying the structure 
function in compensated form, shows the extension of the scaling region. 

The same data have also been fitted by means of the expression proposed in [II] to model the behavior of a generic 
structure function in the entire range of scales (bottom panel, solid line). In our case, the interpolation function for 
5i2(r 3 ) is given by 

1 + B 2 (r 3 /r?) 2 

and describes the superposition of a scaling behavior with coefficient £2(2) at intermediate scales, a large-scale satu- 
ration and a dissipative closure at small scales. Here, the exponent £2(2) is fixed by the direct fit estimated from the 
compensated plot in the top panel of figure 4, the transverse integral scale £t is evaluated according to its definition 
given in table I, while A 2 and B 2 are the only two fitting constants. The ability of equation (5) to correctly reproduce 
the experimental data can be appreciated by looking at the excellent agreement between the solid curve and the open 
symbols in the bottom panel of figure 4. Such an extra fitting procedure, not strictly necessary here, turns out to be 
useful later in the context of the TBL flow. 

Results concerning higher-order statistics of anisotropic fluctuations are reported in figure 5. In the top panel, 
the two transverse observables of order four, namely 51112^3) and 51222(73) are shown, both in their standard and 
compensated forms. A best fit yields for the exponent £2(4) a value of 1.7 ± 0.1, indicated by the solid line. The 
associated error accounts for both the deviation from a pure scaling law and for the slightly different behavior of 
5in2(" , 3) with respect to 51222(73)- The corresponding compensated structure functions of order four are shown in 
the inset of the top panel. Mixed structure functions of order six are displayed in the bottom panel of figure 5. 
In particular, with the configuration of figure 3, three transverse observables can be measured, namely SnimO^), 
5ni222(?"3) and 5122222(^3)- Only the first two, with the best statistical properties, are shown in the figure. Here 
again, minimal differences in the scaling behavior of these quantities are observed. Typically, lower values of the 
exponents are achieved for structure functions with largest weight on the vertical velocity component, i.e. for the 
sixth order S\ 22222 . The exponents increase with increasing weight of u, i.e. for the sixth order moving from 5122222 
through 5ni222 to 5imi2. However, we find that a unique value of £2(6) = 2.05 ±0.15 is able to fit satisfactorily the 
set of statistics of order six. 
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FIG. 4: Homogeneous shear flow. Top: Log-log plot of the mixed structure function of second order, — Si 2 plotted as a function 
of the transverse separation rs. All quantities are normalized with respect to the inner scaling defined by !i T = \Jrjp and 
v ju T . The solid line corresponds to the slope 1.2. The minimal length scale resolved is of the order of the dissipative scale. 
In the inset the same data are shown in compensated form, — S^/ns 1 ' 2 . Bottom: The same data are fitted by means of the 
expression given in equation (5), with A2 = 0.018, B2 = 0.062 and D2 = 3.27, solid line. 



The scaling exponents extracted by the present procedure are consistent with the values given in [4] and in [11], 
namely £ 2 (2) = 1.05 1.22, £ 2 (4) = 1.42 -f- 1.56 and £ 2 (6) = 1.71 + 2.14. 

B. The turbulent boundary layer 

To address the effect of the shear intensity and of different boundary conditions, we consider the more complex 
environment represented by the near-wall region of a fully developed turbulent boundary layer. In this flow config- 
uration two basic difficulties emerge. The first one is associated with the thinness of the region where significant 
changes of the mean gradient occur which poses severe restrictions on the probe dimensions. The second problem 
is related to the relatively large fluctuation level in the lower part of the log-region, which may cause troubles with 
Taylor hypothesis. Concerning the first issue, a boundary layer as thick as possible was realized in a relatively large, 
yet well-controlled, experimental facility, as described in section II. 

Moreover, to ensure a sufficient probe resolution, the measurements were limited to the log-region, sufficiently far 
from the boundary. The second point was instead entirely by-passed by addressing anisotropic observables depending 
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FIG. 5: Homogeneous shear flow. Top panel: The purely anisotropic mixed structure functions of order 4, S1112 (triangles) and 
S1222 (circles). The solid line indicates the slope of 1.7. In the inset the same data are compensated with r 1 ' 7 , same symbols. 
Bottom plot: Two purely anisotropic mixed structure functions of order 6, £111112 (triangles) and 5111222 (circles). The solid 
line indicates the slope of 2.05. In the inset the same data are compensated with r 2 ' 05 , same symbols. 

only on the transverse separation , as already explained for the HS configuration. Also for the boundary layer, the 
data we discuss are presented in dimensionless form, using the friction velocity 

U T = \/T w /p (6) 

as characteristic velocity scale, where t w is the average shear stress at the wall (section II). This corresponds directly 
to the normalization used in subsection III A, since in the near-wall region the total shear stress r is constant in 
the wall-normal direction. The main issue is connected to the assessment of the anisotropic properties at changing 
distance from the wall. 

In figure 6 we summarize the results for the fully anisotropic transverse structure functions already introduced in 
the previous section. The three panels show, from top to bottom, the observables of order 2, 4 and 6, respectively, 
while the different symbols correspond to different distances from the wall, from y + = 350 down to y + = 90. As one 
can see, the scaling properties are not as clear as in the HS case. Independently of the order of the structure function, 
farther from the wall a distinct scaling range emerges. As the wall is approached, the scaling behavior is less evident 
and a tendency towards saturation at large scales is observed. Here, in order to extract quantitative results one needs 
to consider also large scale effects. In particular, we generalized the expression (5) to all orders, for separations much 
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larger than the Kolmogorov scale: 



S ai ... a Jr 3 ) = A n {y+)rl 



1 + Dn(y+)^ 



-6(n) 



(7) 



By comparing the fit with the raw data, it is quite clear that the poor scaling closer to the wall is substantially 
explained in terms of saturation occurring earlier and earlier as the wall is approached. The inset in the top panel 
describes the fitting procedure and highlights the scaling law by removing the effect of the large scale saturation. The 
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FIG. 6: Turbulent boundary layer: log-log plot of the mixed structure functions vs transverse separation r$ at different wall 
distances y + — 350 (triangles), 240 (diamonds), f40 (squares), 90 (circles) and for different orders n = 2, 4, 6. In the top panel 
we show — Si 2 (second order) where the corresponding solid lines represent the fit (7) keeping fixed £2(2) = 1.2 at varying 
A2, D2. The inset shows — S12 in raw form (open triangles) and after compensation with saturation and viscous contributions 
(filled triangles) at y + = 350. The solid line indicates the slope 1.2. In the middle panel we show — S1112 (fourth order), with 
£2(4) = 1.7 at varying A4, D4. In the inset we show the anisotropic flatness — S1112/S12 at corresponding locations. The solid 
line indicates the slope —.65. In the bottom panel we show — S111112 (sixth order) with $2(6) = 2.05 at varying Ae, De- 
inset of the middle panel shows the anisotropic flatness, — Sn22/Si2> to highlight the high degree of intermittency 
showed by anisotropic fluctuations, independently on the distance from the wall. It is important to stress here that 
the good agreement with the data for all distances is obtained by keeping fixed the scaling exponents £2(1) to the 
value obtained in the HS case for all orders n -only prefactors change at changing the wall distance. 

IV. CONCLUDING REMARKS 

In this paper we have performed a systematic analysis of small-scales anisotropic turbulent properties in two different 
experimental set-up, a homogeneous shear and a turbulent boundary layer. We have used two cross-wires probes in 
order to extract the leading anisotropic fluctuations of two-point correlation functions in the homogeneous directions 
without the needing of Taylor hypothesis. We have analyzed structure functions up to order n = 6 finding a good 
agreement of the anisotropic exponents between the two experimental set-up and at changing the distance from the 
wall in the turbulent boundary layer. We could compare the anisotropic properties by changing the normalized shear 
intensity, S 1 *, by a factor 2 and more. 

The cleanest data are obtained for the homogeneous shear flow where a fit of the power law behavior allowed for a 
direct measurement of the anisotropic properties. In the turbulent boundary layer, we had to take into account also 
large scale saturation effects, especially close to the wall, in order to obtain a global fit of the structure functions 
behavior for all value of y + . Our results support the "exponent-only" scenario of universality also in the anisotropic 
sector. In other word, we have been able to fit all experimental data by keeping fixed the scaling properties and 
adjusting only the prefactors. 

The main drawback of all actual experimental set-up consists in the impossibility to exactly disentangling different 
anisotropic sectors among themselves. This implies that sub-leading contributions coming from sector with higher 
j's could spoil the quantitative measurements. In particular, the systematic differences observed here and in other 
studies [4, 11] between the scaling of anisotropic correlation functions of the same order but with different tensorial 

components as S111222 ( r ) an d ^111112 ( r ) ma y weu be due to the effects of sub-leading contributions coming from the 
j > 4 sectors in the SO (3) expansion. 
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